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PREFACE 


This  technical  report  summarizes  work  completed  on  the  project 
for  the  period  May  1972  to  May  1973.  The  following  personnel  of  the 
Amy  Materials  &  Mechanics  Research  Center  contributed  to  the  studies 
in  the  indicated  areas. 

F.  Baratta,  C.  Freese  -  Fracture  Mechanics 
J.  Campo,  6.  Driscoll  -  Mechanical  Tests 

R.  Brockelman  -  Nondestructive  Evaluation 

S.  Der  Boghosian  -  X-ray 

Dr.  H.  Priest  -  Chemical  Vapor  Deposited  Silicon  Nitride 
Dr.  R.  N.  Katz,  A.  King  -  Fractography 

Dr.  E.  M.  Lenoe  -  Statistical  Considerations  &  Report  Preparation 
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1.0  INTRODUCTION 

The  necessity  to  cope  with  temperatures  beyond  the  capabilities  of 
most  metallic  materials  is  the  impetus  to  use  the  so-called  brittle 
ceramics.  Typica,  of  such  applications  are  re-entry  vehicles  (leading 
edge  elements,  nose  caps,  control  surface  panels,  etc.)  and  propulsion 
systems  (rocket  engine  nozzles,  combustion  chambers,  etc.).  Most  recently, 
however,  application  of  ceramics  to  high  temperature  gas  turbine  engine- 
has  provided  additional  stimulus  to  further  develop  design  and  analysis 
methodology  for  brittle  materials. 

.For  examp  1 e ,  Ford  Motor  Company  and  Westinghouse  Electric  Corporation 
have  had  research  programs  underway  for  several  years  on  the  application 
of  ceramics  to  turbine  engines.  In  dune  of  1971.  the  Advanced  Research 
Projects  Agency  of  the  Department  of  Defense  sponsored  a  contract  to 
accelerate  this  wort,  and  to  deaenstrate  that  brittle  materials  can  be 
successfully  utilized  in  demanding  high  temperature  structural  applications 
Two  demonstration  items  were  chosen,  namely  a  vehicular  gas  turbine  and 
portions  of  a  large  turbine  power  generator  with  emphasis  on  the  foraer. 

Thus  far  a  number  of  contractor  reports  have  been  published  on  the 
ARPA-F0RD-WEST1NGH0USE  project(l)*.  In  addition,  many  hundreds  of  hours 
have  been  accumulated  in  hot  engine  test  rigs  on  various  first  and  second 
stage  stators  and  rotors,  shrouds,  nose  cones,  and  combustion  chambers, 
regenerator  seals,  etc.  The  engine  project  is  strongly  oriented  toward 
an  experimental-trial  assessment  of  iterative  designs.  This  also  allows 
proof  testing  of  individual  components,  prior  to  assembly  for  a  required 

^Numbers  in  parenthesis  pertain  to  references 
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200  hour  engine  operation  demonstration  test. 

addition  to  duties  as  contract  unitors  on  this  project,  the 
Army  Materials  4  Mechanics  Research  has  been  conducting  in-house 
research  in  support  of  the  project.  This  report  documents  portions  o 
the  AMMRC  studies  completed  during  the  period  May  1972  to  May  1973.  The 
overall  research  consists  of  theoretical  and  experimental  investrgatron 
of  constitutive  eguations,  fracture  panics,  failure  theories,  thermal 
and  mechanical  fatigue  and  probabilistic  aspects  of  brittle  ,atena 
behavior.  Materials  to  be  investigated  include  several  types  of  si  icon 
nitride  and  silicon  carbide.  Under  ARPA  sponsorship  the  Ford-Westing  ouse 
teara  is  conducting  an  extensive  program  of  mechanical  and  physical 
properties  deteimination  for  the  candidate  ceramic  materials.  However 
much  additional  infection  ought  to  be  developed  and  considerable  effort 
„  required  to  further  refine  design  and  analysis  procedures  for  such 


materi  als. 


2.0  MECHANICAL  PROPERTIES  CHARACTERIZATION 

During  this  project  a  hunter  of  test  methods  are  being  employed 
to  study  response  of  the  ceramic  materials.  Behavior  under  simple 
monotonic  loads  is  being  observed  via  two  types  of  tension  tests  (includ¬ 
ing  hydroburst  of  thin  walled  rings)  and  several  types  of  flexure  loadings 
and  specimen  sizes.  Results  obtained  to  date  have  been  at  ambient  labora- 

tory  conditions  (73°F). 

Tension  experiments  were  completed  on  specimens  of  the  type  shown  in 
-  Figure  1.  These  dogboneconfigurations  are  subjected  to  hydrostatic  pressure 
in  the  apparatus  shown  in  Figure  2.  The  ratio  of  shoulder  to  gage  section 
.  diameter  isselected  to  minimize  the  amount  of  pressure  required  to  fracture 
the  specimens  in  tension.  Furthermore  the  specimen  ends  rfde  on  ”0”-ring 
seals  and  bending  effects  therefore  tend  to  be  reduced  by  this  test  method. 
2.1  TENSION  TESTS  ON  HOT  PRESSED  SILICON  NITRIDE  (NORTON  HS-130) 

Sixteen  specimens  of  the  design  and  method  given  in  AMMRC-TR  69-02 
and  shown  in  Figure  1  were  tested  to  failure.  The  first  specimen  tested 
was  in  the  as-received  condition  such  that  circumferential  (perpendicular 
to  the  load  axis)  score  (grinding  marks  existed  all  along  the  gage  length. 
(See  Figure  3).  Accordingly  the  specimen  failed  at  38.160  psi .  The  next 
specimen  to  be  tested  was  first  carefully  hand-polished  using  diamond  grit 
and  it  failed  at  54,700  psi.  Subsequently  a  simple  polishing  machine  was 
constructed  to  prepare  the  remaining  specimens.  The  resulting  data  is 
tabulated  next. 
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Photomicrograph  of  Typical  Tension  Specimen 


Another  attempt  was  made  to  improve  over  the  as-received  surface 
integrity  by  use  of  chemical  vapor  deoosited  (CVD)  Silicon  Nitride 
coatings  on  the  HPSN.  Two  tension  specimens  were  coated  at  United 
Aircraft  Research  Laboratories.  Hopefully  the  CVD  coating  would  provide 
an  inexpensive  method  of  improving  the  surface  integrity  of  the  HPSN, 
rather  than  costly  polishing  procedures.  Unfortunatley  one  of  the 
specimens  was  destroyed  during  the  CVD  process.  The  other  was 
returned  intact,  however  the  shoulder  region  was  eroded.  Subsequently 
this  specimen  failed  at  33,600  psi.  Inspection  of  the  fracture  surface 
revealed  a  fairly  uniform  deposit  of  silicon  nitride  having  a  nominal 
thickness  of  4  mils.  There  was  an  obvious  Interface  of  diss'milar 
material  and  fracture  had  been  initiated  at  a  flaw  of  nominal  size 
.013  X  .075  inches,  located  just  beneath  the  coating  in  the  mid-portion 
of  the  gage  section.  Obviously  further  development  of  the  CVD  process 
is  required  to  achieve  a  successful  compatible  coating  procedure. 

2.2  FLEXURE  OF  HOT  PRESSED  SILICON  NITRIDE  (Norton  HS-130) 

Beam  bending  experiments  were  conducted  on  HPSN  using  1^  point  )l| 
loading  on  rectangular  specimens  each  3  1/4”  X  0.15  X  0.60  inches.  The 
beams  were  carefully  ground  by  Norton  Company  and  no  extra  surface 
polishing  was  attempted.  The  test  apparatus  is  illustrated  in  Figure  4. 
and  the  resulting  data  is  tabulated  next.  Specimens  were  instrumented 
with  strain  gage  rosettes  both  top  and  bottom  of  the  beam  in  order  to 
measure  poisson's  ratio  and  to  compare  the  tension  and  compression  modulus 
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TABLE  2  SUMMARY 


OF  FLEXURE  TESTS  HPSN 


Modulus 
XI 0-6  psi 

Specimen 


1* 

2* 

3 

4 

5 

6* 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 


standard 


44.44 

44.44 

46.78 

46.18 

45.38 

45.40 

46.47 

45.40 

46.54 

46.00 

45.94 

47.03 

45.82 

46.3 

47.03 

46.36 

46.36 

46.54 

45.82 

44.89 


mean* 
deviation 
variance  % 


46.167 

.568 

1.23% 


*Not  including 


Strength 

ksi 


88 ,889 

88.889 

98.889 
97,222 
74,167 
64,440 
95,277 
82,778 

109,444 
95,556 
100,000 
95,833 
72,222 
85,111 
100,722 
106,667 
101,389 
107,611 
103 ,194 
96,444 

95,443 
10,542 
11  .045% 


1U 


Remarks 

Broke  out  of  g*s. 

II 

Badly  chipped  beam  edges 
Broke  under  roller 

Chipped  edges 

Badly  chipped  edges 


ROUND  BEAM  TESTS 

A  number  of  the  hot  pressed  silicon  nitride  tension  specimens  fractured 
leaving  fairly  large  portions  of  the  straight-sided  gage  section  intact, 
and  it  was  decided  to  perform  flexure  experiments  on  these  fragments. 

Four  point  bend  tests  were  performed  on  eight  specimens  using  the  loading 
span  indicated  i.i  the  sketch.  u 


1.375" 


TABLE  3 

The  experiments  used  rigidly  fixed  upper  and  lower  loading  heads  and  were 

completed  in  a  20,000  lb.  capacity  Instron  using  a  1000  pound  load  cell. 

Specimen  Diameter  Failure  Strength  Remarks 

Load  (lb) 

iq  onn  3?3  59,120  As  machined  finish 


Results  questionable, 
bad  cell  unbalanced 

Break  near  load  point 


Specimen 

Diameter 

Failure 

Load  (lb) 

Strength 

19 

.200 

323 

59,120 

20 

.1975 

518 

98,460 

9 

.198 

722 

136,190 

4 

.1995 

450 

82 ,980 

12 

.1985 

450 

84 ,240 

16 

.200 

358 

65,520 

22 

.198 

422 

79 ,600 

7 

.199 

450 

83,610 

For  specimens  20,  4,  22  and  7  _  m^an  = 

standard  deviation  =  7262 
coefficient  of  variation  =  3.43% 

For  all  specimens  except  #19  niean  =  82,401  psi 

standard  deviation  =  9625 

coefficient  of  variation  =  11.68% 
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23  Reaction  Bonded  Silicon  Nitride  (RBSN)  (Fort  ™erial) 

'  ft  series  of  molds  were  obtained  from  Reliable 
for  injection  adding  of  five  types  of  RBSN  species  using  the  Fort 

.  orocess  The  molds  were  supplied  to  Ford 

Motor  Company  proprietory  process 

.  hnnp.  (See  Figure  1)  and  compact  tension 
and  beams,  rings,  tension  dog  bones  (See 

.  The  ring  specimens  were  unsuccessful  due  to 
plates  were  prepared.  The  ri  g  p 

cracking  associated  with  shrinkage,  the  tension  dogbones  part, a  y 
successful  and -the  bea*s  and  plates  appeared  adeguate.  It  should 
noted  that  the  processing  of  this  material  is  improving  as  ti.e  proc 

and  therefore  the  .suits  reported  herein  are  representative  P  a 
•  i  &<-  »  later  date  additional  specimens  will  be 

generation  material.  At  a  later 

prepared  using  improved  production  techniques. 


♦Reliable  Tool  &  Die  Company, 


Dearborn,  Michigan 


\ 
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Tension  Tests 


The  RBSN  specimen  dimensions  prior  to  testing  are  as  indicated 
elow: 


#  D1 

1  .995 

2  .996 


TABLE  4 


D2 _ d _ 

.9955  -2  to  .188 


(Eccentricity  Tensile 
run-out)  Strength ,psi 

*  7460 


.995 


.184  to  .195  .003  7460 


As  stated  these  specimens  were  injection  molded.  It  appears  that  wnen  cast 
there  was  a  mismatch  in  the  die  which  caused  the  diameters  to  be  egg  shaped. 
Also  run  out  between  diameters  "DJ  and  "d"  was  considerable.  To  correct 
some  of  the  run  out,  diameters  X  were  turned  down  to  a  smaller  size,  built 
up  with  epoxy,  and  .emachined  in  relation  to  diameter  "d".  Because  of 
considerable  mismatch  in  sDecimen  #1,  diameter  "d"  still  appeared  egg 
shaped  and  it  was  impossible  to  obtain  a  true  run  out  reading. 

Due  to  noncircularity  as  well  as  eccentricity,  poor  strengths  were 
achieved.  Examination  of  the  fract.-'&  surfaces  showed  gross  defects  as 
depicted  schematically  in  Figure  5. 


Specimen  #  2 


.0051 


0174 


0249 


0091 


Figure  5.  Schematic  of  Flaws  in  Reaction  Bonded  Silicon  Nitride 

Tension  Specimens. 
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2.4  Flexure  Tests  on  Reaction  Bonded  Silicon  Nitride  (RBSN)  (Ford  Material) 
Flexure  tests  were  also  completed  at  room  temperature  on  RBSN  using 
the  injection  molded  and  nitrided  specimens  supplied  by  Ford.  Both 
rectangular  and  approximately  circular  specimens  were  used.  The  data  is 


listed  next. 


FLEXURE  ON  REACTION  SINTERED  SILICON  NITRIDE 

3.25  X  .15  X  .6  inch  Beams  1/3  point  loading 


TABLE  5 


Specimen 


Modulus  X  10  psi 


Strength,  psi 


1 

13.96 

2 

13.58 

3 

13.80 

4 

14.78 

5 

13.89 

6 

15.65 

7 

15.15 

8 

14.34 

9 

13.58 

10 

14.54 

11 

15.40 

12 

14.58 

13 

14.29 

14 

12.84 

15 

14.19 

16 

14.67 

mean 

17 

14.19 

standard 

deviation 

variance= 


14.32  X  106  =  mean 
.6851=  standard  deviation 


15.556 

12.444 
14,667 
16,000 

13.556 
16,111 
15,333 
16,889 
14,77b 
15,611 
14,778 
.5,111 

16.556 

11 .444 

14,917  psi 
1,474 

iT^HTT  =  .098813  =  9.9% 
14,917 


4,78%  variance 
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round  beams 
.8 


1.375" 


Three  round  beam  specimens  were  also  subjected  to  flexure  tests. 
These  sped  mens  were  obtained  from  tension  specimens  which  were  broken 
unintentionally  prior  to  tensile  testing. The  flexure  loading  just 


described  was  also  used. 

TABLE  6 

Specimen 

1 

2 


3 


Dia 

Failure 

Load  (lb) 

Strength 

.195 

78 

15,400 

.195 

74 

14,610 

.195 

85 

16,790 

All  breaks  occured  in  the  constant  moment  section 


2  *  nondestructive  evaluation 

The  original  .terl.l  billets  as  well  as  final  machined  specimens  have 

,  .  A  T.„  billets  of  hot  pressed  Norton  HS130  silicon 

been  nondestructive^  examined.  Ten  billets 

nitride,  each  six  inches  square  and  nominally  one  and  an  eighth  inch  t  ic 
were  subjected  to  x-ray  radiography  as  well  as  ultrasonic  characterization, 
number  of  the  billets  are  in  the  process  of  being  machined  into  beams  tension 
end  torsion  specimens,  compact  tension  specimens  (figure  6  )  and  notched  be  m 
:l  Of  rupture.,  specimens.  The  available  specimens  have  been  destructively  tested 

as  reported.  In  general  both  the  X-ray  and  ultrasonic  inspection  revea  e 
uniform  quality  material.  So™  small  regions  of  high  density  were  observed 

^  X" Ultrasonic  compressions!  wave  veloicty  measurements  were  made  in  the  plane 

end  perpendicular  to  the  plane  of  the  compact  tension  crach  growth  specimens 

•x-  c x  ,  ciinhtlv  higher  modulus  might 

(See  Figure  6).  The  avenge  velocities  suggest  a  slightly  nig 

exist  along  the  pressing  direction. 

2  5  1  HYDROBURST  SPECIMENS 

In  addition  to  the  attempt  to  produce  NBSN  rings,  three  other  materia  s 
were  considered,  namely  HPSN,  chemical  vapor  deposited  silicon  carbide  and 

silicon  nitride. 

Thirty  of  the  hot  pressed  silicon  nitride  rings  were  obtaine 
Norton  CO.  and  subjected  to  X-ray  radiography,  local  ultrasonic  egression 
and  shear  wave  measurers  and  "C"-scan  observation.  The  x-ray.examination 
served  to  loca.e  several  high  density  regions.  The  "C"-scan  results  appear 
1n  figures  7  and  8.  These  two  illustrations  illustrations  represent  "C  - 
scan  observations  where  the  rings  were  merely  reversed  in  the  test  apparatus. 


( 
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The  nonde.'tructive  evaluation  of  thirty  rings  yield  the  following  defects 
in  twelve  rings 


TABLE  7 

SUMMARY  OF  DEFECTS  IN  RING  SPECIMENS 

Specimen 

B-2 

Nature  of  Defects 

High  density  surface  faw  on  inside  wall 

Y-5 

Two  high  density  areas,  approximately  18 

61 ,  3  and  4 

Outer  wall  high  density  surface  flaws 

W-4 

Hign  density  surface  flaw,  inside  wall 

R-5 

Several  interior  high  density  regions 

0-2 

Outer  wall,  high  density  flaw 

0-3 

Mid  ring  (interior)  high  density  area 

0-1 

Worst  "C"  scan 

Y-4 

Bad  "C"  scan 

6-4 

Bad  "C"  scan 

Y-5 

Bad  "C"  scan,  indications  180°  apart 

Note  that  "C"  scanning  appears  to  be  sensitive  to  low  porosity  regions 
whereas  X-ray  radiography  is  suited  to  determine  high  density  areas.  These 
rings  were  prepared  over  sized  (nominal  dimensions  0.8"  inside  die..  1.35” 
outside  dia.,  X  .625”  high)  in  order  that  the  specimens  might  be  carefully 
machined  to  thin  walled  rings.  During  this  machining  process  a  number  of 
the  detected  flaws  will  be  removed.  However  the  specimens  which  remain 
flawed  will  be  subjected  to  detailed  fractography  in  order  to  establish 
the  significance  of  the  various  flaws. 
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Ultrasonic  through  transmission  measurements  were  completed  on  a  local, 
point-by-point  basis  using  both  compressions!  and  shear  waves.  Each  ring  was 
tested  through  the  wall  thickness  as  well  as  through  the  height  of  the  ring  (which 
corresponds  to  the  pressing  direction  of  HPSN.  The  rings  were  systematically 
scanned  and  the  high  and  low  readings  were  recorded.  Precise  measurements 
of  ring  density  were  also  completed.  Using  these  observations  and  applying 
elementary  formula  of  wave  propagation  theory.  It  was  then  possible  to  estimate 
young's  modulus,  shear  modulus  and  Poisson's  ratio.  The  infection  is  recorded 
1n  Table  8  and  9.  At  a  later  date  the  data  will  be  correlated  with 

destructive  test  results.  At  the  present  time  the  information  can  be  used  to 
estimate  statistical  variability  of  moduli  in  the  Norton  H5-130  material. 


TABLE  8  ULTRASONIC  WAVE  MEASUREMENT  THROUGH  RING  WALLS 


HOT  PRESSED  S^ 


SAMPLE 

Longitude 

Time 

Wall 

Thickness 

Inch 

Density  s 
(gr/cm3' 

Poisson's 

Ratio 

Average 

Young's 

Modulus 

Poisson's 
Ratio 
Aver.  L-H 

Young's 

Modulus 

Y-l 

H* 

L* 

1.27616 

1.24733 

.276 

3.1712 

.2655 

44.84 

46.94 

.2546 

45.84 

Y-2 

H 

L 

1.33864 

1.29589 

.28725 

3.1707 

44.14 

47.10 

.2719 

43.53 

46.45 

Y-3 

H 

L 

1.33864 

1.29588 

.2875 

3.1638 

44.12 

.2547 

45.09 

48.12 

Y-4 

H 

L 

1.34459 

1.30822 

'  .2865 

3.1614 

43.39 

45.84 

.2644 

43.49 

45.94 

Y-5. 

H 

L 

1.37313 

1.24966 

.2865 

3.1731 

41.76 

50.42 

.255 

42.66 

51.51 

0-1 

H 

L 

1.39529 

1 . °4956 

.28585 

3.1725 

40.25 

50.19 

.2807 

38.90 

48.51 

0-° 

H 

L 

1.36063 

1.31347 

.2873 

3.152 

42.49 

45.59 

.2607 

42.91 

46.05 

0-3 

H 

L 

1.28592 

1.26716 

.2752 

3.1621 

43.78 

45.09 

.258 

44.46 

0-4 

H 

L 

1.33487 

1.28218 

.28665 

3.1695 

44.19 

47.89 

.2657 

44.17 

47.87 

0-5 

H 

L 

1.34191 

1.27108 

.287 

3.1617 

44.78 

48.89 

.2781 

43.54 

47.88 

B-l 

H 

L 

1.26710 

1.24649 

.27425 

3.1617 

44.78 

46.27 

.2781 

43.54 

44.99 

B-2 

H 

L 

1.26307 

1.25002 

.274 

3.166 

45.04 

45.99 

.2753 

44.08 

45.01 

B-3 

H 

L 

1.34718 

1.29102 

.2745 

3.1679 

39.76 

43.30 

.2533 

40.75 

44.37 

B-4 

H 

L 

1.26498 

1.25362 

.27425 

3.1715 

45.07 

45.89 

.2788 

43.75 

44.55 

B-5 

H 

L 

1.35035 

1.27945 

.2875 

3.1702 

- 

43.44 

48.39 

.2599 

43.95 

48.96 

♦NOTE:  The  rings  were  examined  on  a  point-by-point  basis  and  the  high  (H)  and  low  (L) 
values  were  recorded.  23 


n_l  H  1.36430 

K  L  1.26922 

.286k5 

3.1778  .2655 

42.26 

48.83 

.2564 

43.05 

49.75 

r_2  H  1.25542 
"  '  L  1.22914- 

.26915 

3.1521 

43.80 

45.69 

.2642 

43.92 

45.82 

45.44 

45.88 

„  3  H  1.23516 

K  L  1.22916 

.269 

3.1623 

45.34 

45.79 

.2645 

R  4  H  1.34405 

K*  L  1.30708 

.2875 

3.1746 

43.91 

46.43 

.2635 

44.10 

46.63 

n  5  H  1.34402 

K  L  1.29289 

.28595 

3.1767 

43.47 

46.98 

.256 

44.32 

47.89 

u.l  H  1.36168 

.28685 

3.1778 

42.63 

49.93 

.2538 

43.65 

51.13 

L  1.25821 

u_ 2  H  1.33420 

W  L  1.29615 

.2868 

3.172 

44.31 

46.95 

.2627 

44.57 

47.23 

44.54 

44.71 

u  3  H  1.26367 

W  L  1.26127 

.27415 

3.193 

45.43 

45.61 

.2746 

W-4  H  1.34221 

.2869 

3.1786 

43.91 

47.18 

.2755 

42.95 

46.15 

L  1.29485 

42.27 

48.66 

W-5  H  1.37400 

.2862 

3.1799 

41.71 

48.02 

.259 

L  1.28058 

G-l  H  1.35075 

.287 

3.1725 

43.30 

47.55 

.2594 

43.85 

48.15 

L  1.28892 

G-2  H  1.36075 

.286 

3.1776 

42.44 

48.52 

.2788 

41.20 

47.10 

L  1.27263 

G-3  H  1.36900 

.28635 

3.1606 

41.80 

M  .1*1 

.2735 

41.08 

46.89 

L  1.28137 

G-4  H  1.37902 

.28675 

3.1824 

41.60 

49.39 

.2801 

40.26 

47.80 

L  1.26557 

6-5.  H  1.26017 

.2732 

3.1566 

44.85 

45.42 

.258 

45.55 

46.13 

L  1.25224 
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TABUE  9  ULTRASONIC  WAVE  MEASUREMENTS  THROUGH  RING  HEIGHT 
HOT  PRESSED  Si 3N4  RINGS 

“2"“*  W)  W*  Velocity 


SAMPLE  Time 


Poisson's 

Ratio 


Y-l 

H  2.65710 

L  2.65113 

4.62573 

0.5798 

3.1712 

46.74 

46.80 

436416 

437398 

250685 

250685 

i 

.2538 

.2554 

Y-2 

H  2.64583 

L  2.64274 

4,72440 

0.5794 

3.1707 

45.37 

45.40 

437972 

438484 

245280 

245280 

.2715 

.2723 

1 

.2539 

.2554 

Y-3 

H  2.66029 

L  2.65469 

4:63189 

0.5797 

3.1638 

46.49 

46.55 

435817 

436736 

250308 

250308 

I 

.2599 

.2688 

Y-4 

H  2.69118 

L  2.65673 

4.72475 

0.5799 

3.1614 

44.90 

45.21 

430963 

436552 

245473 

245473 

Y-5 

H  2.66045 

4.63124 

0.5800 

3.1731 

46.69 

46.78 

436016 

437574 

250473 

250473 

.2537 

.2563 

L  2.65098 

• 

0-1 

H  2.64314 

L  2.63200 

4.77664 

0.5798 

3.1725 

44.75 

44.84 

438721 

440578 

242765 

242765 

.2793 

.2820 

.2538 

.2675 

Or  2 

H  2.71305 

L  2.66115 

4.72322 

0.5798 

3.152 

44.56 

45.05 

427416 

435751 

245510 

245510 

.2524 

.2635 

0-3 

H  2.69664 

4.68558 

0.5794 

3.1621 

45.31 

45.71 

429720-' 

436435 

247312 

247312 

2.65515 

0-4 

H  2. *7089 

L  2.65085 

4.71086 

0.5796 

3.1695 

45.53 

45.53 

437294 

437294 

246070 

246070 

.2683 

.2683 

0-5 

H  2.65684 

4.76079 

0.5796 

3.1722 

44.82 

44.89 

436308 

437801 

243489 

243489 

.2738 

.2761 

2.64778 

R-l 

H  2.64806 

L  2.64090 

4.62009 

0.5797 

3.1778 

47.00 

47.07 

437830 

439017 

250947 

250947 

.2554 

.2573 

.2583 

.2701 

R-2 

H  2.70202 

L  2.65643 

4.73339 

9.5796 

3.1521 

44.50 

44.92 

429012 

436375 

244898 

244898 

.2619 

.2671 

R-3 

H  2.6643 

L  2.64666 

4.69493 

0.5795 

3.1623 

45.49 

45.68 

434664 

437910 

246862 

246862 

R-4 

H  2.63531 

4.64174 

0.5800 

3.1746 

46.81 

46,91 

440176 

441819 

249906 

249906 

.2622 

.2648 

L  2.62551 

* 

.2551 

.2568 

{ 

R-5 

H  2.65650 

4.63275 

0.5801 

3.1767 

46.78 

46.84 

436740 

437777 

250434 

250434 

L  2.65021 

25 
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Longitude  Shear  Thickness 


SAMPLE 

Time 

Time 

Inches 

B- 1  H 

L 

2.63741 

2.62639 

4.75695 

0.5799 

B-2  H 

L 

2.65294 

2.64576 

4.75748 

0.5797 

B-3  H 

L 

2.66040 

2.65530 

4.62381 

0.5800 

B-4  H 

L 

2.63574 

2.62125 

4.74581 

0.5799 

B-5  H 

L 

2.64299 

2.62327 

4.62263 

0.5796 

W-l  H 

L 

2.65334 

2.64445 

4.61107 

0.5798 

W-2  H 

L 

2.62989 

2.61715 

4.62447 

0.5781 

W-3  H 

L 

2.64495 

2.64231 

4.74204 

0.5799 

W-  H 

L 

2.64279 

2.63829 

4.74337 

0.5800 

W-5  H 

L 

2.62492 

2.60990 

4.58950 

0.5757 

G-l  H 

L 

2.64739 

2.64273 

4.64061 

0.5792 

G-2  H 

L 

2.63970 

2.62281 

4.75084 

0.5799 

G-3  H 

L 

2-.  7031 4 
2.65717 

4.79999 

0.5793 

G-4  H 

L 

2.64776 

2.63734 

4.73100 

0.5798 

G-5  H 

L 

2.68196 

2.65346 

4.67106 

0.5797 

Density  Young's  longitude  Shnar  Poissons 

(gr/cnr)  Modulus  Velocity  Velocity  Ratio 


3.1617 

44.93 

45.03 

439750 

441595 

243812 

243812 

.2781 

.2807 

3.166 

44.82 

44.89 

437025 

438211 

243700 

243700 

.2744 

.2761 

3.1679 

46.72 

46.77 

436025 

436862 

250875 

250875 

.2526 

.2540 

3.1715 

45.25 

45.37 

440028 

442461 

244384 

244384 

.2770 

.2805 

3.1702 

46.88 

47.08 

438594 

441891 

250766 . 
250766 

.2572 

.2625 

3.1778 

47.09 

47.18 

437034 

438503 

251482 

251482 

.2525 

.2550 

3.172 

46.77 

46.90 

439638 

441778 

250018 

250018 

.2610 

.2644 

3.193 

45.53 

45.55 

438496 

438934 

244578 

244578 

.2742 

.2749 

3.1786 

45.34 

45.38 

438930 

439679 

24452 

244552 

.2749 

.2760 

3.1799 

47.06 

47.21 

438642 
441 1 66 

250877 

250877 

.2569 

.2610 

3.1725 

46.55 

46.59 

437563 

438335 

249622 

249622 

.2588 

.2600 

3.1776 

45.23 

45.37 

439368 

442197 

244125 

244125 

.8767 

.2808 

3.1606 

43.67 

44.06 

428613 

436028 

241376 

241376 

.2678 

.2791 

3.1824 

44.78 

44.81 

437955 

439685 

242543 

242543 

.2788 

.2813 

3.1566 

45.62 

45.90 

432296 

436939 

248209 

248209 

.2541 

.2618 

2.5.2  Reaction  Bonded  jn»™n  Nitride  Flexure  Specimens 

The  RBSN  flexure  speciemsn  previously  described  were  also  subjected 

to  nondestructive  evaluation.  A  total  of  fifty-six  beams  were  provided  by 
Ford  for  both  room  and  elevated  temperature  experiments.  Those  specimens 
were  porous  (~75%  theoretical  density)  and  had  obvious  surface  flaws 
on  one  face  of  the  beam.  Typical  "C"  scan  observations  are  shown  in 
Figure  9  for  both  sides  of  the  samples.  Apparently  the  lines  correspond 
to  the  furnace  rack  supports  used  during  the  nitriding  cycle.These  are  visually 
appanrfon  only  one  face  and  accordingly  this  poorer  surface  was  used  on 
the  compression  side  during  flexure. 
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2.5.3  ULTRASONIC  VELOCITY  MEASUREMENTS 

ON 

REACTION  SINTERED  SILICON  NITRIDE  STATORS 

Two  ceramic  stators  were  subjected  to  ultrasonic  examination. 
Compression  waves  were  introduced  through  both  si  das  of  the  rim  of  the 
stator  wheel  near  the  blade  root  area.  Velocity  transit  times  were 
recorded  for  each  location.  At  a  later  date  these  areas  will  be  subjected 
to  mechanical  testing  and  localized  density  measurements.  The  flexural 
modulus  and  flexure  strength  will  be  determined  using  miniature  beams 
machined  from  the  area  of  each  blade  root.  Similar  test  will  be 
conducted  on  beams  prepared  from  the  individual  stator  blades.  While 
the  strength  measurements  are  not  expected  to  be  representative  of  actual 
mechanical  properties  of  the  reaction  sintered  silicon  nitride,  the  data 
is  expected  to  be  representative  of  spatial  variation  of  moduli  and 
strengths  throughout  these  production  parts.  Note  the  data  presented 
In  the  following  tables  suggest  high  uniformity  in  modulus  throughout 

the  parts. 


TABLE  10  STATOR  NUMBER  1 


SITION 

THICK. 

TIME 

TIME 

Velocity 

Velocity 

(IN) 

TL 

TS 

VL  in/sec 

VS  in/ sec 

A-l 

.102 

.682705 

1.15093 

298,811 

177,248 

A- 2 

.10 

.679185 

1.13017 

294,471 

176,965 

A- 3 

.099 

.668972 

1.11757 

295,977 

177,170 

A- 4 

.099 

.661205 

1.10182 

299,433 

179,703 

A- 5 

;}§05 

.661740 

1.10764 

302,234 

180,564 

A- 6 

.658878 

1.10542 

305,064 

1 81 ,831 

A- 7 

•  TO 

.668772 

1.12205 

299,056 

178,245 

A- 8 

•  1005 

.674270 

1.13099 

298,100 

1 77 ,720 

A- 9 

.101 

.677813 

1.13287 

298,017 

1 78,308 

A- 10 

.10105 

.68880 

1.13225 

293,375 

178,494 

A-l  1 

.102 

.692256 

1.15626 

294,689 

176 ,431 

A- 12 

.1032 

•  .696000 

1.15955 

296,552 

1 78 ,000 

A- 13 

•  1034 

.706389 

1.17512 

292,757 

175,982 

A-l  4 

.1002 

.705612 

1  .15071 

284 ,009 

174,153 

A- 15 

.102 

.703661 

1  .16426 

289,912 

175,219 

A- 16 

.102 

.704699 

1.16945 

289,485 

174,441 

A-l  7 

.102 

.707977 

1.16647 

288,145 

174,887 

A- 18 

.104 

.703791 

1  .17261 

295,542 

177,382 

A- 19 

.105 

.686914 

1  .16992 

305,715 

179,449 

A- 20 

.105 

.689339 

1.16745 

304,640 

179,879 

A-21 

.105 

.688489 

1.15258 

305,016 

182,200 

A- 22 

.1035 

.694473 

1.16624 

298,068 

177,493 

-23 

.1025 

.677188 

1.14715 

302,722 

178,704 

A- 24 

.102 

.664549 

1 .120021 

306,975 

1 82  ,’|40 

A- 25 

.101 

.679029 

1  -15154 

297,484 

175,417 

B-l 

.103 

.704320 

1.18809 

292,481 

173,388 

B-2 

.  103 

.714875 

1.20585 

288,162 

170 ,834 

B-3 

.103 

.702683 

1.17587 

293,162 

175 ,189 

B-4 

.103 

.682249 

1.15621 

301  ,943 

178,168 

B-5 

.1035 

.687424 

1.16393 

301 ,124 

1 77 ,846 

B-6 

.1025 

.691121 

1.16003 

296,626 

176 ,720 

B-7 

.102 

.697800 

1.9323 

292,347 

1 70,965 

B-8 

.1015 

.694337 

1.7737 

292,365 

172,418 

B-9 

.1013 

.679722 

1.15788 

298,063 

174,975 

B-10 

.0995 

.683448 

1.16884 

291 ,171 

170,254 

B-ll 

.099 

.659644 

1.11760 

300,162 

177 ,165 

B-l  2 

.098 

.669816 

1.14853 

294,111 

171 ,524 

B-l  3 

.0976 

.664173 

1.15253 

293,899 

169,367 

B-l  4 

.0968 

.666322 

1.162011 

290,550 

166 ,608 

B-l  5 

.097 

.658538 

1.11898 

294,592 

173,372 

B-l  6 

B-l  7 

B-l  8 

.0968 

.632487 

1.10686 

306,093 

174,909 

B-l  9 

.097 

.659145 

1.11795 

294,321 

173 ,539 

B-20 

.098 

.674405 

1.1 3608 

290,627 

172,523 

B-  21 

.099 

.677398 

1.16048 

292,295 

170,619 

-22 

.099 

.681805 

1.16258 

290 ,406 

170,311 

d-23 

.10 

.682806 

1.16923 

292,909 

171 ,053 

B-24 

.1012 

.692123 

1.19066 

292,434 

169,990 

B-25 

.1016 

.6/4874 

1.16322 

30 

301  ,093 

174,688 

TABLE  11  STATOR  NUMBER  2 


PARTS 

THICK. 

(IN) 

LONG 

SHEAR 

Vl  in/sec 

Vs  in/sec 

Microseconds 

C-l 

.257 

1.56574 

2.68215 

328.279 

191 ,637 

Notch 

C-2 

.249 

1.56734 

2.75043 

317,735 

181 ,063 

c  3 

.252 

1.55490 

2.79278 

324,137 

180,465 

High  Att. 

C-4 

.2475 

1.53038 

323,449 

C-5 

.250 

1.52177 

2.62011 

328,565 

190,832 

C-6 

'.247 

1.52132 

2.60980 

324,718 

189,287 

C-7 

.251 

1 ,55483 

2.66285 

322,864 

188,520 

C-8 

.251 

1 ,53662 

2.64056 

326,691 

190,111 

C-9 

.2475 

1 ,53849 

2.65446 

321 ,744 

186,479 

C-10 

.247 

1 ,54184 

2.65211 

320,396 

186,267 

Hiqh  Att. 

C-ll 

.252 

1  ,55776 

2.70057 

323,541 

186,627 

C-l  2 

.2413 

1  ,53291 

2.62459 

314,826 

183,876 

D-l 

.1975 

1  ,20977 

2.06419 

326,508 

191,358 

D-2 

.197 

1,21678 

2.08911 

323,305 

188,597 

D-3 

.202 

1 ,22440 

2.12636 

329,958 

189,996 

D-4 

.205 

1 ,27909 

2.16658 

320,540 

189,238 

D-5 

.204 

1 ,24884 

2.15082 

326,703 

189,695 

D-6 

.204 

1 ,25778 

2.14086 

324,381 

190,578 

D-7 

.1965 

1,21106 

2.08497 

324,509 

188,492 

D-8 

.198 

1.21971 

2.08298 

324,667 

.  190,112 

D-9 

.200 

1.22976 

2.08813 

325,267 

191,559 

D-l  0 

.2015 

1,21832 

2.13440 

330,783 

188,812 

D-l  1 

.1965 

1.21622 

2.07586 

323,132 

189,319 

D-l  2 

.203 

1.23816 

2.11275 

327,906 

192,167 

\ 
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3.0  PROBABILITY  BASED  DESIGN  &  ANALYSIS 

The  basic  idea  of  probability  based  design  and  analysis  is  to  attempt 
to  design  a  specified  reliability  directly  into  a  component,  rather  than  to  use 
a  conventional  safety  factor  approach  which  is  based  merely  on  average  values 
of  strengths  and  mean  load  conditions,  Benefits  of  the  approach  include  the 
ability  to  assess  the  sensitivity  c:  the  design  to  variability  in  material 
properties,  mechanical  and  thermal  loads,  geometric  tolerances,  etc.  Further¬ 
more  the  probabilistic  approach  provides  an  estimate  of  failure  rates  of  the 
components  of  the  design,  under  expected  service  conditions. 

There  are  a  variety  of  methods  to  incorporate  statistical  considerations 
into  stress  analysis.  In  this  report,  we  begin  first  with  determination  of 
appropriate  statistical  models  to  represent  the  strength  behavior  of  the 
ceramics.  In  particular  the  room  temperature  tensile  strength  of  hot 
pressed  silicon  nitride  is  usedin  determining  suitability  of  various  simple 
models. 

WEIBULL  DISTRIBUTION 

Consider  a  material  which,  under  a  uniform  stress  condition,  has  the 
probability  of  failure  of  P£  per  unit  volume.  Then  the  overall  probability  of 
survival  for  volume  V  of  the  material  is  given  by 
Ps  =  0-Pf)V 

for  the  simultaneous  survival  of  V  volume  units.  Then 
ln  S  =  V  In  (1-Pf) 

Weibull  statistics  (3)  defines  the  risk  of  rupture  as  R  =  -ln  S  and  in  terms 
of  an  infinitesimal  element,  gives 
dR  =  -ln(l-Pf)  dV 


\ 
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where  the  term  -ln(l-Pf)  Is  assumed  to  be  a  positive  function ,  * fa)  of  the 

tensile  stress  O'  ,  alone.  Substitution  gives  dR  -  r.((T)  dV 

R  -  ft dv 


From  this  definition  of  R  we  arrive  at 
Ps  =  exp 


lltion  OT  R  VIC  ai  I  M 

[- 1 


Several  for  -  for  n  (O')  were  proposed  by  Weibull,  e.g. 

(  iZL  \rrn 
n  (CD  =  ^  Co ' 

•  w  cCT-^-r 


and  in  these  cases 
P 


“t-  fv  izr«i 


Ps  = 


«*  [-rr-^rav] 


Another  form  was  proposed(4)  to  account  for  a  maximunvalue  of  stress, 
above  which  the  probability  of  survival  vanishes 

cr  -co“Uj  O’  >crs 


where 


/infer)  -  o 


EFFECT  OF  SIZE  AND  STRESS  DISTRIBUTION 

Using  the  simplest  form  of<y><$  ^  we  comPare  tw0  types  of  sPecimens 
under  uniform  stress  states  at  the  same  value  of  failure  probability  then 

the  ratio  of  stresses  for  the  two  volumes  is 

23  .  ( )k)'u 
O*  ' v  V/ 

Strength  is  therefore  shown  dependent  on  volume  and  decreases  as  the  volume 
increases.  Furthermore  the  effect  becomes  less  important  as  the  value  of  the 
"Weibull"  modulus  m  increases.  The  equation  can  be  generalized  in  the  form 

where  K  is  a  normalized  geometry  integral  and  CTis  a  reference 
stress  state.  For  the  same  type  of  test  specinens,  the  size  effect  will  be  as 
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theoretically  defined.  Thus  by  performing  experiments  on  various  sites  of 
specimens  and  also  by  using  different  types  of  loadings,  it  is  possible 
to  compare  and  assess  the  applicability  of  such  representations  for 
observed  material  strengths.  First,  however,  we  begin  by  apply  ng  several 
Of  the  more  co«on  methods  to  fit  VleibulTmodels  to  experimental  data. 


V 


3.1  EVALUATION  OF  WEI  BULL  PARAMETERS 

1,1.1  THEORETICAL  CURVE  FITTING  METHODS 

METHOD  OF  MOMENTS  -  Three  Parameter  Wei  bull 
The  method  of  moments  is  based  on  choosing  the  Wei  bull  parameters 

i,e*  ovij  ^ 


so  that  the  first  three  moments  of  the  distribution  correspond  to  those 
of  the  observed  data.  The  first  step  is  to  calculate  these  3  moments  froi 

the  data.  Let  us  label  these  as 
Oyyn  =  Mean  strength 
G?~  =  Square  of  the  standard  deviation 
n  “Third  moment  about  the  mean 


Furthermore,  let  Vj  =  Number  of  datapoints  and. 


0^  *  The  strength  of  the  "n"  specimen 


Tben  the  first  moment 

(T  m\  - 


Second  moment 

qa-' 


Third  moment 


1 

)  <Tnn 

(1) 

N 

J_ 

(2) 

N 

l  ■ 

n\:\ 

J_ 

1 1 

^“(cr-<r»p 

(3) 
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Consideration  of  the  theoretical  distribution  and  taking  the  corresponding 
moments  of  all  points  on  the  probability  curve  about  the  smallest  point 
or,  leads  to  the  following  theoretical  relations 

or,  5k..  r( \+-h) 


(4) 


where  KV  is  the  appropriate  constant  K  and  Volume  V  for  the  particular 
type  of  test  specimen  and  is  the  Gamma  Function.  Furthermore 


a*--  [r(l+^l-  r’^|+i)] 

(KV)™ 


(5) 


r  =  £  [r( i+^)-3rO^)r(i^)4'2P,(t^)]  (6) 

(KV)^ 


These  equations  could  be  solved  by  trial  procedures  for  /mnt  0? }  cru  • 
However  this  would  be  cumbersome.  Instead  by  dividing  equation  (6)  by  (5) 
raised  to  3/2  power ^we  define 

A *3  “  (tf ) ^ 

v‘^y  (which  can  be  computed  from  the  data) 


nfl+js:) 

-3  PI 

r( 

1 

1  +  o.  r  (\  ■+■  ^ 

[r( 

,1  +  ^. 

0  + 

m> 

(7) 


This  new  "variable"  A3  can  therefore  be  calculated  as  a 
function  of  Wei  bull  parameter  m.  Either  a  graphical  procedure  or  a 
computer  search  technique  can  be  used  therefore  to  estimate  the  value 
of  m  which  corresponds  to  the  calculated  A^. 
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Now  by  rearranging  equation  (5)  we  can  also  define 

Al=  [r(i*a-'’^*>J  » 

end  now  A2  may  bellculated  using  the  value  on.  estimated  via  A, 

Then  since  <X  (M™  1*  a  known  quantity,  OJ  is  calculate 

equation  (8).  Finally  let  A,-  r’0+rvn) 

then  we  can  reyrite  (4)  a 

rr-  Go  r\\ 


<Tu  -  ^  - 


(Kv)’a 


These  equations  have  been  programmed  and  Iterative  procedures  adopted 
to  apply  this  technique.  Examples  of  the  method  appear  next. 


TABLE  12 


EVALUATION  OF  PARAMETERS  BY  METHOD  OF  MOMENTS 
HOT  PRESSED  SILICON  NITRIDE 


;rengthCTn,  Ksi 

CM 

C 

b 

(JV 

85.41 

7294.87 

623,054.68 

78.77 

6204.71 

488,745.23 

76.77 

5893.63 

452,454.20 

73.54 

5408.13 

397,714.00 

69.820 

4874.83 

340,360.80 

65.46 

42851 .01 

280,496.9 

65.39 

4275.85 

279,597.97 

63.62 

4047.5 

257,502.23 

63.02 

3971.52 

250,285.21 

62.88 

3953.89 

248,620.88 

61.82 

3821.71 

236,258.26 

54.70 

2992.09 

163,667.32 

52.67 

2774.13 

146,113.37 

(y  =  i  y  (Tm 

=  67,221 

Mean  strength 

<r«.)  =  W1* 


Square  of  Standard  Deviation 


T"T" 

\. 


TABLE  13 


itandard  Deviation 

85.41 

18.189 

78.77 

11.549 

75.77 

9.549 

73.54 

6.319 

69.82 

2.599 

65.46 

1.761 

65.39 

1.831 

63.62 

3.601 

63.02 

4.201 

62.88 

4.341 

61.82 

5.401 

54.7 

13.151 

52.67 

14.551 

330.84 

133.3794 


91.1834 

39.9298 

6.7548 

3.101121 

3.352561 

12.9672 

17.6484 

18.84428 

29.1708 

172.9488 

221.7316 


a 


' 2 . 


(1081.852163)  =  83.2194 
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V 


Third  Moment 


M 


i  V”  o"  3  _  3  ^  fa  •+  2 

-  N  L  N  Z- 


M  = 


2  5  £» X 10 


M 


A,  =  -2.  =  .337 

W  /rn  ^2-6 


At  ^  '  l35- 

or  =  74Si 


Ki  —  '88T 

on = 


and  supposedly  the  probability  of  failure 


l  -  evf 


L 


A 


2-Z51 
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3.1.2  GRAPHICAL  PROCEDURES 

The  starting  three  parameter  eou^  ca^ be  written  as 

=  e*p  ^  } 

I  ~ri 

or  taking  logarithms  /  la/  \ 

n  ctai-Ufical  distribution  function  will  be  linear 

Now  theoretically  the  statist! ca  ,  _  \ 

if  the  data  is  plotted  on  coordinates  ^^(/-|/ersus  ™  P 
and  the  slope  of  the  linear  curve  will  yield  the  Weibull  modulus  .. 


An  example  of  evaluation  of  the  Keibull  parameters  is  presented  next 

TABLE  14 


Rank 

Probability  of 
Failure  n 

N+l 

1  nl 

Strength 

psi 

Two  Parameter 
Distribution 

Q«\  <r 

Three  Parameter 

(o r-os ; 

1 

.07143 

-2.602 

52,670 

10.8718 

8.9189 

2 

.1429 

-1 .870 

54,700 

10.9096 

9.1593 

3 

.2143 

-1.422 

61 ,820 

11.03198 

9. 7185 

4 

.2857 

-1 .089 

62,880 

IT. 04898 

9.7803 

5 

.357 

-0.817 

63,020 

11.0512 

9.7882 

6 

. .4286 

-0.580 

63,620 

11.0607 

9.8213 

7 

.500 

-0.366 

65,390 

11.0881 

9.9130 

8 

.5714 

-0.165 

65,460 

11.0892 

9.9165 

9 

.6429 

0.0292 

69,820 

11.1537 

10.1114 

10 

.7143 

0.225 

73,540 

11.2056 

10.2521 

,1 

.7857 

0.432 

76,770 

11.2486 

10  3600 

12 

.857 

0.666 

78,770 

11.2743 

10.4215 

13 

.9287 

0.970 

85,410 

11.3552 

10.6019 

No  -  that  the  value  =  45, 1 23  was  obtained  previously  by  the  method  of  moments. 
Otherwise  a  trial  procedure  would  have  been  employed  to  arrive  at  va.ue  of  0^  • 


The  data  for  the  two  and  three  parameter  distribution  are  shown 
in  Figure  10  Comparison  of  the  graphical  procedure  and  the  method 
of  moment  models  is  made  in  Figure  H  .  It  is  apparent  that  there  is 
a  certain  amount  of  flexibility  in  fitting  the  rather  limited  data; 
and  several  models  could  justifiably  be  chosen  to  represent  the 

measurements. 

NORMALIZED  STRESS  METHOD 

In  this  method  (6)  the  equations  suggested  by  Weibull  are  rewritten 
with  the  mean' strength  substituted  for  and  0"u  ^  • 

Then  Ps  can  be  written  as  ^ 

Ps  e«p-L^r<^] 

where  and  O'  Is  the  failure  stress 

associated  with  a  particular  probability  of  survival  P$. 

The  procedure  for  evaluation  of  a  set  of  d£ta  by  this  method  is  to 
normalize  the  observations  by  dividing  by  the  mean  strength  and  rank 
these  in  order  of  increasing  size.  As  examples  of  the  method,  the 
normalized  strength,  ,  and  corresponding  probability  of  survival 

values  are  tabulated  for  tension  and  flexure  data  on  HPSN  and  flexure 
data  for  RBSN.  The  values  are  plotted  in  Figuresl2,  13.  Note  that 
in  the  region  ft  near  1.0  a  small  uncertainty  is  probability  results 
in  a  wide  band  of  m  values.  If  the  data  tend  to  a  vertical  line  then 
a  single  two  parameter  Weibull  distribution  describes  the  data.  This 
is  apparently  not  the  case  for  the  test  results  reported  here.  (See 

Figure  12). 
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Rank  R 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 
11 
12 
13 


TABLE  15  TENSION  STRENGTH ,  HPSN 


Strength 

Probability 

R-.J 

52,670 

.077 

.784 

.948 

54,700 

.154 

.815 

.873 

61 ,820 

.231 

.922 

.798 

62 ,080 

.308 

.938 

.724 

63,020 

.385 

,942 

.650 

63,620 

.462 

.947 

.575 

65 ,390 

.538 

.973 

.500 

65,460 

.615 

.976 

.425 

69 ,820 

.692 

1.04 

.35 

73,540 

.770 

1.097 

.277 

76,770 

.846 

1.142 

.200 

78,770 

.923 

1.17 

.127 

•15,410 

1.00 

1.27 

.052 

N  -  13  mean  =  67,200 

Refer  to  Figure  11,  and  disregarding  three  lowest  values,  m  =  8.94 


TABLE  16  FLEXURE  STRENGTH,  HPSN 

i 

ank  R 

Strength 

79 

U-M 

1 

0 

72,722 

.764 

.95 

74,467 

.778 

.9023 

C. 

3 

82,778 

.870 

.845 

W 

A 

85,111 

.895 

.787 

5 

95,277 

1.0 

.73 

V 

6 

95,556 

1.001 

.673 

V 

7 

8 

q 

95,833 

96,444 

1.012 

.615 

.557 

97,222 

1.02 

.500 

10 

98,889 

1.037 

.442 

11 

100,000 

1.05 

.385 

12 

100,722 

1.125 

.327 

13 

101 ,389 

1.062 

.270 

14 

103,194 

1.085 

.214 

15 

106,667 

1.117 

.155 

16 

107 ,611 

1.13  • 

.098 

17 

109,444 

1.15 

.022 

N  =  17 

Disgarding 

mean  = 

values  near  B  =  1, 

95,443 

then  from  Figure  12, 

m  =  10. f 

47' 
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TABLE  17  FLEXURE  TESTS  ON  RBSN 


Strength 

Rank 

.  _ 

1  N+-4 

& 

11,444 

1 

.9513 

.77 

12,444 

2 

.882 

.835 

13,556 

3 

.812 

.91 

14,4667 

4 

.743 

.984 

14,778 

5 

.673 

.99 

14,778 

6 

.613 

.99 

15,111 

7 

.535 

1.013 

1533 

8 

.465 

1.03 

15,556 

9 

.398 

1.042 

15,611 

10 

.327 

1.047 

16,000 

11 

.257 

1.073 

16,111 

12 

.188 

1.082 

16,556 

13 

.118 

1.11 

16,889 

14 

.048 

1.132 

mean  =  14,917 

and  from  Figure  12,  m  =  14.3 


CUMULATIVE  SURVIVAL  PROBABILITY 


There  are  a  number  of  alternate  approaches  to  describe  the  results 
if  the  data  are  curved  about  either  side  of  a  vertical  distribution. 
First  the  three  paramter  distribution  might  better  represent  the  data. 
Alternatively  a  combination  of  two  parameter  distributions  might  better 
fit  the  measurements.  In  the  case  of  the  three  parameter  distribution, 
/0  is  defined  as  T/""  •  The  tension  strengths  for  HPSN 

are  evaluated  this  way  in  the  following  tabulation  and  in  Figure 
The  value  of  CTU determined  via  the  method  of  moments  is  chosen  to  treat 

the  data,  e.g.’Q^=  45,198 


TABLE  18 


Rank 

0"; 

a 

i  -  £ 

1 

52,670 

.3393 

.948 

2 

54,700 

.4315 

.873 

3 

61 ,820 

.7548 

.798 

4 

62,880 

.8029 

.724 

5 

63,020 

.80924 

.65 

6 

63,620 

.8345 

.575 

7 

65,390 

.9169 

.5 

8 

69 ,460 

.920 

.425 

9 

69 ,820 

1.118 

.35 

10 

73,540 

1.287 

.277 

11 

76,770 

1.434 

.2 

12 

78,770 

1.524 

.127 

13 

85,410 

1.826 

.052 

N  =  13 

and  the  range  in  m  from  Figure  13  is  1.2  to  3.8 
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From  the  preceding  discussion,  it  appears  that  for  the  rather  small 
sample  sizes  used  here,  the  normalized  stress  method  is  not  useful 
for  values  nearfrl  .  Other,  and  simpler  model  fitting  techniques 
ere  advisable.  It  is  interesting,  however  to  compare  the  results 
obtained  from  the  different  techniques  just  discussed. 


TABLE  19  SUMMARY  OF  WEIBULL  PARAMETERS 
OBTAINED  BY  VARIOUS  TECHNIQUES 

Two  Parameter  Model 


Three  Parameter  Model 
2.55 

2.08 


Method  of  Moments 

51,163 

8.9  7491 

45,198 

Graphical 

46,931 

6.85  6009 

45,198 

Normalized 

8.94 

Distribution 

From  the  earlier  discussion  we  have  that  Mear^ Fail u^ 

Strength  GTtri  -  c|(J“ 

°u 


Variance 


C*=  fVckc>)  -CX"' 

*  cru 


2.34 


and  the  risk  of  rupture 

R  -  T 


52 


V 


r t  anH  0™, 

Suppssedly  the  three  distribution  parameters  m, 
independent  of  size  and  manner  of  loading  and  in  general 

R-KVt^r 

where  the  load  factors  have  been  presented  for  simple  mechanical 

tests  by  a  number  of  authors  (6,7). 

For  a  two  parameter  model,  the  mean  failure  strength 
rr  CTo  rU4^) 

°7m.-  “ 


U  are 


and 


C^~  0~o 

^  =  (KV)^ 


[rfi.ii-rtiviJ] 


and  under  simple  tension  K  the  loading  factor  becomes  1.0.  .Figure 
14  shows  the  estimated  mean  strength  volume  dependency  for  HPSN  in 
comparison  to  test  data.  Considering  data  variability  the  models 
appear  reasonably  consistent.  When  additional  data  is  available 
the  analytical  expressions  will  be  further  refined  and  surface,  as 
opposed  to  volume  dependency  will  also  be  surveyed. 


AMMRC  Tension  Data 
Ford/Westinghouse  Data 


Weibull  Fit  to  AMMRC  Flexure  Data 


Method  of  Moment 

Graphical  Fit  to 
AMMRC  Tension  Data 


LO  Volume,  inches' 


R  Tensile  Strength  versus  Volume  Hot  Pressed  Silicon 

Nitride,  Norton  HS  130 


3.2  Choice  of  the  Probability  Model 

Suppose  the  test  specimens  could  be  assumed  representative  of  the 

production  parts  in  all  respects,  including  the  surface  machining.  Then 
our  test  data  can  be  assumed  to  represent  statistical  variation  of  the 
material  properties  and  there  remains  the  question  of  choosing  the  appropriate 
statistical  model  to  fit  the  data. 

Quite  often,  for  instance,  material  strength  properties  are  taken  to 
have  normally  distributed  properties  (8).  Given  the  fact  that  test  data  is 
ordinarily  too  limited  in  quantity,  and  this  does  not  allow  a  precise 
theoretical  discrimination  between  various  analytical  models,  the  actual 
statistical  model  is  often  selected  on  the  basis  of  convenience. 

Reliability  theories  have  been  fairly  extensively  developed  for 
normal,  log-normal,  two  and  *hree  parameter  tleibull  and  other  statistical 
distributions.  Each  distribution  has  unique  qualities  and  can  be  applied 

in  various  ways. 

Thus  unless  there  are  firm  theoretical  basis  for  choosing  one  particular 
distribution,  the  selection  is  made  after  applying  simple  goodness  of  fit 
criteria  and  ease  of  use  of  the  data.  The  normal  distribution,  for  instance 
is  easily  applied  to  simple  design  formula.  Haugen  (9)  has  discussed  two 
-  pararreter,  statistical  algebra  wherein  elementary  reliability  calculations  are 
stated  in  terms  of  mean  values  and  standard  durations  of  the  variables. 
However,  before  discussing  rrethods  of  applying  other  distributions,  let  us 
consider  several  of  the  properties  of  the  Weibull  distribution. 


;  ) 
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RELATIONSHIP  BETWEEN  PROBABILITY  BASED  STRENGTH 

AND 

FACTOR  OF  SAFETY 

Suppose  the  material  strength  can  be  represented  as  a  two  parameter 
Weibull  distribution  then  the  Probability  of  Survival,  Ps  is  given  by 

P,  -  \  -  ff  =  e*p  [-  ^ 


or 


fs  =  [-  r 


^vn 

with  0"  =  applied  stress 

=  mean  strength 

Then  let  SF 


=  conventional  safety  factor 

<r 


and  in  effect  this  previous  equaticn  is  a  relation  between  probability 

safety  factor  and  Weibull  modulus,  m. 

The  safety  factor,  SF,  can  also  be  expressed  as 

<r 


sf  - 


[e*  (. 


_L  \  “I  Y*ir\ 


rj- 1 


in  which  for  high  reliability  levels,  the  approximations 

for  Pp  <  10 


Pf+®2+  -  -  P 

Pp=|-PS  , 

Then  SF  ^ 

and  since  T  (  l  +  ^  “  *  for  most  values  of 

fs  *  (& 


m 


yrri 


irm 
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Now  if  a  structure  is  composed  of  00  subvolumes,  each  with  a  volume 
equal  to  that  of  the  test  specimen,  or  alternatively  using  tne 
appropriate  Wei  bull  parameters  then  the  safety  factor  for  the  array 
of  subvolumes  under  uniform  stress  &  is  \ 

r  _  (T/wi  _  P  ( 1  +  ^ 

v  ”  rr  “  n  t  \/  \  V**  ^  ft  ' 


or 


u  -i-p’  -  '”(£) 


rrc 


The  extreme  safety  factor  s  also  of  interest.  This  is  defined 

as  the  ratio  of  the  probable  least  strength  of  the  material 

to  the  probable  design  stress  CTp 

Op-  _  0  im\ 

Using  analytical  methods,  we  can  predict  the  most  likely  least  value 
in  a  set  of  /Y\  specimens  (10).  In  particular  the  probability  of  fY\ 
specimens  simultaneously  surviving  the  stress  0"  is  given  by 

P,  -  LPsCtf] 


<r\ 


>'T» 


Probability  density  is  found  by  d  fferentiating  once  e.g. 

4l9U  d(3 


ft,  =  JE- 

CPm 

The  most  probable  least  value  corresponds  to  the  maximum  of  the  probability 
density  which  is  obtained  by  another  differentiation  and  equating  to  zero. 

™s  y1elds  u  0  ( r, " 

1  ^  Aft  ~  40 
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Using  the  two  parameter  Weibull  expression  leasts  to  the,  most  likely 


From  these  elementary  considerations,  it  is  apparent  that  the  Weibull 
distribution  permits  a  variety  of  useful  general izatiors  and  a  formal 
approach  particularly  suited  for  rapid  assessment  of  reliability  in 
structures  under  simple  stress  conditions  and  where  the  variability 
in  loading  is  not  large,  so  that  predominant  statistical  variation  is  due 
to  strength  as  opposed  to  applied  stresses. 


3.3  APPROACHES  FOR  PROBABILITY  BASED  STRESS  ANALYSIS 

Suppose  we  ere  dealing  with  a  material  which  does  not  exhibit  volume 
dependent  strength.  Furthermore  assume  that  no  variation  in  the  critical 
mechanical  and  thermal  loadings  is  expected.  Let  the  probability  density 


distribution  of  strength  be 


crC  s') 


Let  the  probability  density  distribution  function  associated  with  the  applied 
stresses  throughout  the  structure  and  resulting  from  the  design  load  be 


(TU) 

Then  the  reliability  i.e.  the  conditional  numerical  probability  that  the 
component  will  perform  without  failure  may  be  formally  stated  as 

Reliability  -  R  -  f  CT(S)] dis¬ 

integrated  forms  of  this  mobility  equation  are  available  for  normal ,  lognormal 
gamma,  extreme  value  and  other  common  statistical  distributions.  Refering 
to  the  normal  distribution,  the  reliability  estimate  is  essentially  simplified 
to  computation  of  the  mean  values  and  standard  deviations  of  the  stress  and 
strength  distributions.  Then  with  the  applied  stress  and  strength  probability 


X 


=  mean  values 
=  standard  deviations 
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The  reliability  is  simply  calculated  from _ 

R-  ' 

Cs  6 

where  C\  —  S  —  ^ 

6  =  s-±_ 

Kececioglu  (It  has  also  described  a  transform  method  to  evaluate  the 
general  reliability  expression. 

R  =  j"°°Qfrf'(s)dsJcru')^ 

~oO  ^ 

Introducing  the  notation 


6  = 


F  = 


/  ar( s)  d  s 

jL- 

r<*> 

J  <j(A.)d^ 


Then 


R  =  f  (SdF 

O 

where  dF=  Cf(A  Ja 


and  unreliability  is  expressed  as 
,c>0_  ,~.o C> 


q.  =  r  [  r  d'cAAijtft'&is 

-oO  S 

by  letting  G(s)~  f  cr(s>ds 


.0*5 

FM  = 


Then 


where 


Q  ■=.  J  F 
d(,  -  o-f?)olS 
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Alternatively  the  reliability  expression  can  be  written 

=  /"t  f 


-oo 


-cO 


and  letting 


c(>\  -  f  creels 

1  '  -OO 

GCs)  ro«Js 


-  Oo 


Then 


-  f  F 


where 


and 


da  -  <rCs>as 


where 


q  .  f°T  fo-^dsjcrw^ 

Jo6 

S(s  )=f  CT(fs)ds 


»-  oto 


>4, 


F  _  /  crMcU- 


-  oo 


so  that 


with 


Q  =  /fid  F 

0 

dp  =  Q'(sl)cU- 
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Thus  in  applying  this  method,  a  principal  decision  concerns  selection 
of  the  appropriate  statistical  forms.  This  is  obviously  based  on  the 
strength  data  and  knowledge  of  the  stress  distribution.  Note  that  the 
technique  does  not  incorporate  volume  dependency;  nor  variation  of  material 
moduli,  loadings,  or  structural  geometry. 

determination  of  stress  distributions 
monte  carlo  technique 

It  is  well  known  that  stress  is  generally  a  function  of  many  variables,  i.e 

cr  -  o\  K  Ti 

where  L  through  k  represent  loading,  temperature,  geometry,  material  properties 
time,  stress  concentration  factors,  etc.  The  Konte  Carlo  Technique  is  aimed 
at  determining  the  influence  of  these  variables  on  the  stress  distribution 


function. 

A  brief  summary  of  the  technique  described  by  Kececioglu (12 )  is  as  follows 
1.  Divide  each  distribution  of  the  variables  into  a  number  of  intervals, 


\ 
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2.  Calculate  the  centroids  ;  L  J  j  J 

for  all  the  intervals  in  each  distribution. 

3.  Determine  the  prohabiiity  of  occurence  of  the  centroid  of  each 

interval ,  ft  j  ,  Which  is  the  percenter  of  the  A  distribution  curve 
% 

contained  in  the  J*tVv  interval 

4  Enter  all  the  pairs  of  hunters,  Ci;  ,  and  P-  ,  into  a  digital 
computer,  along  with  a  randOT  number  generating  program  and  a  program 
for  associating  particular  digits  in  the  random  numbers  with  a  partrcu  ar 

pair  of  C:-  and  Pj,  values  of  that  variable. 

5  Generate  a  random  number  and  identify  a  complete  set  of  Ctj  va  ues 
for  all  variables  associated  with  it,  together  with  their  associated 

probabilities  P^i  . 

6.  Calculate  the  stress  value  fr»  each  set  of  the  randomly  selected 

variables. 

7.  For  each  value  of  stress,  determine  the  product  of  the  interva 

probabilities  associated  with  the  C^j  values  used. 

8.  Determine  a  suitable  number  of  intervals  for  the  calculat 

values  and  group  these  stresses  and  probabilities  into  suitable  interva  s. 
calculate  the  centroid  stresses  for  each  interval i  these  are  abscissas 

of  the  <TU)  distribution. 

10.  Determine  the  probability  of  occurance  associated  with  the  central 
stresses  by  s»»ing  up  the  probabilities  in  the  interval.  These  are  the 

ordinates. 

„  These  resulting  distributions  now  are  fitted  by  statistical  regression 
analysis  (normal,  log  nomal,  Weibullland  the  distribution  with  the  highest 

correlation  factor  is  chosen. 
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The  discussion  so  far  has  used  a  maximum  principal  strength 
No  consideration  has  been  made  of  the  appropriate  failure  theories 
for  combined  stresses.  However,  regardless  of  the  failure  theory 
applied,  the  formal  logic  outline  here  is  relevant. 

It  is  apparent  that  a  number  of  formal  statistical  treatments 
can  be  applied  to  estimate  the  reliability  of  a  subcomponent.  Choice 
of  the  appropriate  method  is  intimately  related  to  results  obtained  under 
simulated  service  conditions.  Data  currently  being  obtained  on  the  ARPA- 
FORD-WESTINGrfOUSE  Ceramic  Turbine  project  will  provide  the  information 
to  discriminate  between  these  treatments. 


4.0  FRACTURE  MECHANICS  CONSIDERATIONS 
Fatigue  &  Fracture 

Fracture  mechanics  approaches  emphasize  visible,  geometrically  defined 
crack  and  the  resistance  of  a  given  material  to  its  sloe  or  sudden  growth 
under  load.  Sudden  crack  growth  phenomena  is  based  on  a  well  defined  theory, 
particularly  for  brittle  materials,  whereas  slow  growth  relies  rather  heavily 
on  procedures  of  experimental  materials  evaluation.  Which  is  to  say  that  the 
theoretical  approach  based  on  the  instability  condition  uniquely  defines  the 
start  of  fast  propagation  whereas  slow  crack  growth  dependence  on  stress 
intensity  range'is  often  a  purely  empirical  correlation,  and  subject  to 
significant  statistical  variation.  Nonetheless  fracture  mechanics  procedures 
have  provided  an  engineering  approach  to  deal  with  cracked  or  flawed  structure 
and  allow  generally  conservative  estimates  to  be  made  of  service 

The  usual  application  of  linear  elastic  fracture  mechanics  involves 
two  major  steps.  First  of  all  a  conventional  stress  analysis  is  conducted 
of  the  structure  under  consideration,  ignoring  the  presence  of  defects. 
Usually  this  detailed  stress  information  is  then  coupled  with  simplified 
stress  intensity  calculations  using  mathematical  models  which  resemble  the 
actual  situation  with  respect  to  geometric  boundaries,  stress  distribution 
and  defect  characteristics.  Successful  utilization  of  fracture  mechanics 
theories  is  also  dependent  on  adequate  information  concerning  materials 
properties,  existing  defects,  as  well  as  knowledge  of  stress  conditions. 

The  following  Aerials  properties  are  ordinarily  necessary  for  various 
fracture  mechanics  considerations: 
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-  plane  strain  fracture  toughness  . 

-  conventional  tensile  yield  strength 

-  Young's  modulus 

-  Poisson's  ratio 


crack  growth  rate  (increase  in  crack  length  per 
fatigue  cycle)  as  a  function  of  the  stress 
intensity  and  stress  history 


The  tasic  toughness  parameter  K,c  is  essential  in  all  considerations, 
from  a  simple  comparison  of  materials  to  a  c»plex  calculation  of  allowable 

defect  sizes. 

It  should  be  recognized  that  Kic  is  influenced  by  many  factors.  For 
instance  can  be  measured  via  a  variety  of  test  methods  and  this  can 
introduce  additional  scatter  in  the  observations.  Hence  it  is  essential 
that  the  designer  be  certain  that  the  test  material  condition,  notch 
orientation  and  test  procedures  used  for  fracture  toughness  measurements 

be  fully  representative  of  the  component  of  interest. 

A  number  of  examples  are  available  in  the  literature  to  estimate 
the  importance  of  numerous  variability  factors  for  metals  and  ceramscs. 
These  include  the  recent  evaluation  of  fracture  mechanics  data  for  the 
Mil  aircraft  (13).  Swanson  and  Gross  (14)  for  example,  have  reported 
variation  in  the  order  of  8*  for  K,c  values  for  alumina  oxide,  whereas 
variation  of  twice  this  magnitude  are  not  uncommon  for  certain  ceramics. 

Fracture  toughness  of  materials  is  measured  by  a  variety  of  test 
rethods.  Three  of  the  co«n  types  of  experiments  include  the  compact 
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tension  (CT)  surface  flaw  (SF)  and  double  cantilever  beam  (DCB)  types 
of  specimens.  Figure  15  illustrates  these  configurations. 

Room  temperature  crack  growth  tests  were  performed  at  AMMRC  on 
hot  pressed  and  reaction  sintered  silicon  nitride  using  a  form  of  compact 
tension  specimen.  Only  initial  fracture  loads  could  be  measured  as  the 
crack  rapidly  propagated  across  the  specimen.  The  contoured  double  cantilever 
form  was  also  tried  at  room  temperature  with  similar  results.  Both  configurations 
will  be  used  at  elevated  temperature  to  determine  whether  the  technique  is 

suitable. 

Thus  far  slow  crack  growth  data  has  been  reported  by  Evans  S  Wiederhorn  (15). 
They  employed  an  edge  notched  plate,  under  condition  of  edge  twist  and  for 
constant  displacement  rate  conditions.  The  specimen  is  shown  schematically 

in  Figure  16. 


DISC  TESTING  FOR  FRACTURE  TOUGHNESS  DETERMINATION 
We  have  seen  that  fracture  toughness  of  materials, 
inherent  property,  can  be  experimentally  determined  by  a 


considered  as  an 
variety  of  methods 


such  as  notched  tension  specimens,  notched  flexure  beams,  edge  notched  plates. 
Bursting  of  a  notched  or  otherwise  pre-flawed  disc  can  also  be  used  to  measure 
crack  behavior  under  biaxial  stress  conditions.  Actual  determination  of  this 
fracture  toughness  is  valid  provided  that  brittle  fracture  occurs  and  also 
provided  that  an  accurate  analysis  exists  for  the  experimental  method. 

In  their  early  paper,  Winne  and  Wundt  (16)  assumed  that  the  Bowie 
solution  (17,  18)  for  an  infinite  plate  containing  a  notched  hole,  and 
the  stress  field  in  a  rotating,  bored  and  notched  disc  would  have  equivalent 
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'  CONFIGURATION  AND  GEOMETRY  OF 

CONFIGURATION  AND  GEOMETRY  OF  COMPACT  FLAW  SPECIMEN 


,re  15  Typical  Crack  Growth  Specimens  Compact  Tension,  Double 
r  ‘  ^  pantilpupr  Rparn  anrl  Surfacfi  Flaw. 


stress  distribution  around  the  holes  and  cracks,  provided  the  overall  crack 
length,  L  satisfied  the  condition  ^  r\ 

0,  +  it  <  a 

„here  0,  -  diameter  of  disc  bore  and  02  =  outer  diameter  of  disc.  This 
criteria  was  chosen  by  analogy  to  the  case  of  an  infinite  flawed  plate  as 
contrasted  to  a  finite  width  center  flawed  plate.  Analytical  results  for  | 

the  finite  sized  disc  have  been  completed  by  Freese  and  this  permits 
examination  of  the  validity  of  the  stated  criteria,  as  well  as  reassessment. of 
the  use  of  discs  to  measure  fracture  toughness.  The  results  of  the  stress 
intensity  computation  are  shown  in  Figure  17.  Using  the  same  formal 
procedure  outlined  by  Winne  t,  Wundt,  a  notch  geometry  can  be  chosen  to 
Induce  brittle  fracture  and  the  corresponding  value  of  K,c  can  be  determined. 

The  analytical  results  indicate  the  marked  effects  of  geometry  on  the 
stress  intensity,  and  it  is  obvious  that  the  fracture  behavior  of  the  disc  is  depen¬ 
dent  on  notch  geometry.  Thusitis  not  possible  to  evaluate  or  compare  materials 
characteristics  on  the  basis  of  failure  behavior  of  notched  discs  of  arbitrary 

geometry. 

In  order  to  observe  fracture  toughness  of  the  material,  it  is  necessary 
to  cause  failure  via  brittle,  crack  induced  fracture  at  stresses  considerably 
below  the  tensile  strength  of  the  specimen.  Given  the  fracture  toughness 
of  the  material  and  the  solution  for  stress  Intensity,  we  can  equate  this 

to  maximum  stresses  induced  in  a  thin  rotating  disc. 

The  procedure  would  allow  comparison  of  fracture  toughness  results  under 
loading  rates,  biaxial ity  and  temperatures  similar  to  the  turbine  engine 
environment. 


V 
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Figure  17.  Stress  Intensities  For  Notched  Disc. 


Flaw  Size  Reduction 

There  is  obviously  potential  substantial  benefit  associated  with  studies 
of  fracture  initiating  flaws  in  the  materials.  It  is  reasonable  to  assume 
that  low  strengths  will  be  enhanced  by  some  factor  proportional  to  the 
square  root  of  the  flaw  size.  Furthermore,  if  the  source  of  the  fracture 
initiating  flaws  can  be  determined,  then  appropriate  fabrication  controls 
can  be  instituted  to  effect  the  desired  strength  enhancement.  In  order 
to  apply  fracture  mechanics  considerations,  it  must  be  established  that 
the  crack  extension  conditions  for  flaws  that  control  strength  are  similar 
to  gross,  geometrically  defined  cracks  on  which  stress  intensity  computations 
are  based.  Barring  this  possibility,  the  next  obvious  requirement  is  to 
at  least  investigate  how  well  elementary  fracture  mechanics  predictions 
apply  to  these  parti cualr  ceramics. 

These  computations  must  be  based  on  detailed  fractography  of  the  specimens. 


4.1  FRACTOGRAPHY 

Fractographic  rnalysis  can  be  used  for  a  variety  of  purposes.  Two 
conmon  applications  are  a)  to  determine  the  location  and  size  of  fracture 
origin  in  order  to  provide  data  for  fracture  mechanics  analysis,  and  b) 
to  elucidate  the  nature  of  the  defects  causing  a  fracture  to  initiate, 
relate  this  flaw  to  the  material  process,  then  use  this  information  to 
improve  the  production  method  such  that  critical  flaws  are  eliminated, 

diminished  in  size  or  frequency  of  occurance. 

During  this  program,  fractographic  analysis  was  completed  in  both 
of  these  areas.  As  mechanical  test  specimens  having  clearly  defined 
fracture  origins  become  available,  they  will  be  examined  in  the  Scanning 
Electron  microscope  (SEM)  with  respect  to  geometry,  microstructure  (and 
microtopography)  as  well  as  for  chemical  compositions  of  various  regions 
or  phases  of  interest  using  the  non-dispersive  X-ray  (NDX)  unit  which  is 
attached  to  the  SEM. 

Of  the  fractures  examined  to  date,  fractures  have  originated  both 
internally  and  on  the  surface  of  the  specimens.  When  fractures  initiate 
within  the  central  portion  of  the  bars,  the  initiating  flaw  is  the  intrinsic 
strength  limiting  flaw,  and  not  a  possible  machining  artifact.  Thus  the 
hot  pressed  silicon  nitride  (HPSH)  specimens  with  interior  fracture  origins 
were  subjected  to  special  scrutiny. 

GEOMETRIC  CONSIDERATIONS 

The  HPSN  tension  and  flexure  specimens  described  in  Tables  and 
were  carefully  examined  and  typical  flaw  geometries  were  documented.  A 
typical  tensile  fracture  face  for  interior  failure  origin  appears  in  Figure  *8. 


\ 


73 


IhN  xzi 


t  1740X  Magnification 


while  schematics  of  the  flaw  geodes  observed  to  date  are  shown  in 
Figures  19,r20,  21.  Similar  results  for  the  flexure  specimens  are  shown 
1n  Figure  22.  In  idealizing  these  flaws,  the  gross  shape  was  determined 
at  fairly  low  magnifications  since  observation  of  the  details  using  the 
scanning  ele  ;tror  microscope  tended  to  obscure  the  shape  of  the 


COMPOSITION  ANb  SOURCE  OF  FLAWS 

Examination  of  the  fracture  surfaces  of  available  specimens  indicates 
that  there  are  flaws,  inclusions  or  second  phase  particles  larger  than  the 
initiating  flaw.  Therefore  a  hey  problem  to  be  addressed  in  the  fractographic 
portion  of  this  program  is  identification  of  the  nature  of  the  fracture 
initiating  flaws  and  how  these  differ  from  the  non-fracture  initiating  materia 
defects.  Resolution  of  this  matter  will  undoubtedly  provide  insight  into 
procedures  to  improve  strength  and  the  distribution  of  strength. 

All  flaws  have  been  examined  and  characterized  chemically  by  non- 
dispersive  X-ray  techniques  (NDX)  and  the  composition  of  the  flaw  is 
compared  to  the  average  material  composition.  In  practice  the  NDX  unit 
on  our  scanning  electron  microsocope  is  sensitive  to  elements  with  atomic 
numbers  higher  than  Mg's.  Therefore  N2,  02  or  Mg  is  not  detected  in  our 
routine  studies.  The  equipment  should  detect  Mg.  but  it  is  in  fact  lost 
in  nosie  at  the  end  of  the  scale.  Note  that  the  vertical  axis  of  the  NDX 
traces  in  Figures  23  and  25  have  been  considerably  expanded,  resulting  in  a 
noisy  pattern.  The  instrument  actually  reads  out  electronically  the  Kev 
associated  with  a  given  peak.  Thus  the  ambiguity  suggested  in  these  figures 

is  not.  in  fact  present. 


( 
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Specimen  Number 
#  17  L 


Typical  Defect  Sizes 


Specimen 

#  4  R 

#  4  L 

#  11  R 

#  12  L 

#  12  R 
.  #16  L 

#  16  R 

#  20  R 

#  20  L 


.0007,  .0009,  .0010,  .0012 

.0005,  . 

0007  . 

0009,  . 

0013 

.0024, 

.0013 

.0019, 

.  0016, 

.0010 

.0009, 

.001 

.  0011, 

.0014 

.0005, 

.0006, 

.  0007, 

.  0011, 

.0012 

.0009, 

.  0011, 

.0013 

,.0013 

.0011, 

.0013, 

.0015, 

.0017 

.0007, 

.0011, 

.  0016, 

.  0018, 

.0019 

.0006,  .0007,  .0008,  .0010 
Typical  Defects  in  Tension  Specimen,  Na  4,11,12,16  and  20 


Figure  2L 
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,  fnnr  interior  fracture  origins  indicates  that: 

A  detailed  examination  of  four  interior 

•  in  regions  of  incomplete  densification  with 

a)  Fractures  initiate  in  regions 

, Coo  c-imire  23  and  a  and  b). 

Urge  amounts  of  interconnected  porosity  t 

b)  inclusions  or  second  phases  do  not  apparently  initiate 

Such  areas  of  high  porosity  have  been  observed  to  be  the  initiation  s,  es 

f„r  fracture  in  other  ceramics.  In  the  hot  pressed  silicon  nitride  exa^n 
for  fracture  c  B,  ancl  („  some 

here,  there  is  a  tendency  for  elements  such  as  .  .  '  fracture 

instances  Fe  to  preferentially  segregate  to  the  areas  in  wh, 

originates.  In  second  phase  particles  which  were  not  fracture  mitia  ion 

s1t  s  A1  Ho,  Ca,  and  Fe  have  been  observed  preferentially.  However. 

S  1  ’  A  ’  .  u  ore-'ature  to  attribute  the 

based  on  the  limited  available  data  so  far 

fracture  origin  to  any  specific 
high  porosity-low  density  regions  of  fracture 

impurity  or  conbination  of  impurities. 
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